Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 22 June 2011 



(MN MfcX style file v2.2) 



The large-scale magnetic field and poleward mass accretion 
of the classical T Tauri star TW Hya 

^ i J.-F. Donati 1 *, S.G. Gregory 2 , S.H.P. Alencar 3 , J. Bouvier 4 , G. Hussain 5 , 
O ;M. Skelly 1 , C. Dougados 4 , M.M. Jardine 6 , F. Menard 4 , M.M. Romanova 7 
Y.C. Unruh 8 & the MaPP collaboration 

^ 1 1 IRAP-UMR 5277, CNRS & Univ. de Toulouse, 14 Av. E. Belin, F-31400 Toulouse, France 
1—5 | 2 California Institute of Technology, MC 249-17, Pasadena, CA 91125, USA 

, 3 Departamento de Fisica - ICEx - UFMG, Av. Antonio Carlos, 6627, 30270-901 Belo Horizonte, MG, Brazil 

4 IPAG-UMR 5274, CNRS & Univ. J. Fourier, 414 rue de la Piscine, F-38041 Grenoble, France 

5 ESO, Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany 

6 School of Physics and Astronomy, Univ. of St Andrews, St Andrews, Scotland KY16 9SS, UK 

7 Department of Astronomy, Cornell University, Ithaca, NY 14853-6801, USA 

8 Astrophysics Group, Blackett Laboratory, Imperial College London, SW7 2AZ, UK 



(N 



00 



43 : 

Oh 1 2011 June, MNRAS in press 

6 ' 



> 

(N 

o 



ABSTRACT 

We report here results of spectropolarimetric observations of the ~8 Myr classical 
T Tauri star (cTTS) TW Hya carried out with ESPaDOnS at the Canada-France- 
Hawaii Telescope (CFHT) in the framework of the 'Magnetic Protostars and Planets' 
(MaPP) programme, and obtained at 2 different epochs (2008 March and 2010 March). 
Obvious Zeeman signatures are detected at all times, both in photospheric lines and 
in accretion-powered emission lines. Significant intrinsic variability and moderate ro- 
tational modulation is observed in both photospheric and accretion proxies. 

Using tomographic imaging, we reconstruct maps of the large-scale field, of the 
photospheric brightness and of the accretion-powered emission at the surface of 
TW Hya at both epochs. We find that the magnetic topology is mostly poloidal and 
axisymmetric with respect to the rotation axis of the star, and that the octupolar com- 
ponent of the large-scale field (2.5-2.8 kG at the pole) largely dominates the dipolar 
component. This large-scale field topology is characteristic of partly-convective stars, 
supporting the conclusion (from evolutionary models) that TW Hya already hosts 
a radiative core. We also show that TW Hya features a high-latitude photospheric 
cool spot overlapping with the main magnetic pole (and producing the observed ra- 
dial velocity fluctuations); this is also where accretion concentrates most of the time, 
although accretion at lower latitudes is found to occur episodically. 

We propose that the relatively rapid rotation of TW Hya (with respect to AA Tau- 
like cTTSs) directly reflects the weakness of the large-scale dipole, no longer capable 
of magnetically disrupting the accretion disc up to the corotation radius (at which 
the Keplerian period equals the stellar rotation period). We therefore conclude that 
TW Hya is in a phase of rapid spin-up as its large-scale dipole field progressively 
vanishes. 

Key words: stars: magnetic fields - stars: formation - stars: imaging - stars: rotation 
- stars: individual: TW Hya - techniques: polarimetric 



1 INTRODUCTION 

Magnetic fields play a significant role throughout the life 
of stars, all the way from the cradle to the grave (e.g., 
iDonati fc LandstreetJ |2009| . for a review) . The impact of 
* E-mail: donati@ast.obs-mip.fr magnetic fields is strongest during the very first stages of 
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their life, when protostars and their surrounding plane- 
tary systems form out of the collapse of giant molecular 
clouds; in particular, fields are presumably efficient at slow- 
ing down the cloud collapse, at inhibiting the subsequent 
fragmentation and at dissipating the cloud angular momen- 
tum through magnetic braking and the magnetised outflows 
and colli mated jets generat ed as part of the formation pro- 
cess fe.g jAndre et al.ll2009l ). In the latter phases, when low- 
mass protostars are still actively collecting mass from their 
surrounding accretion disc, the large-scale magnetic fields 
that the protostars generate (through dynamo processes) 
are strong enough to disrupt the inner disc regions, funnel 
the disc material to the stellar surface and drastically brake 
the rotation of the protostar (see, e.g.. lBouvier et al1 l2007a, 
for a review). 

Magnetic fields at the surfaces of cTTSs were first re- 
ported indirectly, through the detection of various spectral 
proxies, usually continuum or line emission throughout the 
whole electromagnetic spectrum, from X-rays to radio wave- 
lengths. It is only about 2 decades ago that fields of cTTSs 
were first detected directly, i.e., through the Zeeman broad- 
ening of spectral lines, and found to reach typical m ag- 
netic intensities of several kG (e.g., Ijohns- Krull 2003, for 
an overview). However, the actual large-scale topologies of 
these fields, controlling in particular how the fields couple 
the protostars to their discs, redirect the accreted disc ma- 
terial and subsequently brake the stellar rotation, remained 
unclear until recently. Thanks to the advent of sensitive 
high-resolution spectropolarimeters dedicated to the study 
of stellar magnetic fields, it is now possible to monitor po- 
larised Zeeman signatures of cTTSs in various spectral lines 
(e.g., those probing accreting and non-accreting regions at 
the stellar surface) and use such time-series to reconstruct 
the parent large-scale magnetic topologies of cTTSs (follow- 
ing the principles of tomographic imaging) ; this option offers 
in particular a more quantitative way of studying magneto- 
spheric accretion processes of cTTSs. 

This is the precisely the main goal of the interna- 
tional MaPP (Magnetic Protostars and Planets) project, 
designed to investigate (through a first survey of a dozen 
targets) how the large-scale magnetic topologies of cTTSs 
vary with stellar par ameters such as mass, age, rotation and 
accretion rates (e.g.. iDonati et al.|[20icl . l2011ah . MaPP also 
aims at providing an improved theoretical description of 
the observational results, using both analytical modelling 
and numerical simula tions (see, e.g., iGregorv et al.l l2010l: 
iRomanova et al]|201l[ ) - the ultimate ambition being to ob- 
tain a deeper understanding of magnetic field generation 
and magnetospheric accretion processes in cTTSs and of 
their impact on the formation of low-mass stars. Results 
up to now demonstrate in particular that the large-scale 
magnetic fields of cTTSs closely resemble those of low-mass 
main sequence dwarfs (once both samples are matched ac- 
cording to their internal structure), with high-mass fully- 
convective cTTSs (with masses ranging from about 0.5 to 
1.0 Mq typically, for ages of a few Myr) hosting strong ax- 
isymmetric dipoles, while partly-convective cTTSs (of either 
higher masses and/or ages) exhibit a more complex field and 
a much weaker (potentially non-axisymmetric) dipole com- 
ponent. 

The new study presented in this paper focusses on the 
evolved cTTS TW Hya, whose relative proximity to the 



Earth and intense magnetic fields l|Yang et al.ll2005l , l2007h 
make it an obvious target o f investigation for MaPP. With 
an age of about 8 Myr (e.g.. iTorres et al.ll2008T ) and a mass 
of about 0.8 Mq (see Sec. [2]), TW Hya can be viewed as 
an evolved version of 2 other prototypical cTTSs, namely 
AA Tau and BP Tau, and can thus serve as a landmark 
to test new generation models of stellar formation including 
the effect of magnetic fields. After summarising the main 
characteristics of TW Hya, (Sec. [5J, we describe the new 
spectropolarimetric observations we collected (Sec. [3} and 
outline the associated temporal variability and rotational 
modulation (Sec.|3|. We then detail the modelling of these 
data with our magnetic imaging code (Sec. [5}, briefly com- 
pare our results with those published already and conclude 
with a updated summary of what MaPP results tell us about 
how magnetic fields impact the formation of Sun-like stars 
(Sec.©. 



2 TW HYA 

TW Hya (K6V) is one of the low-mass cTTSs located near- 
est to the Earth (at a distance of 57 ± 7 pc); it belongs to 
a loose association of young stars after which it was named 
(the TW Hya association), whose age is estimated to be 
~8 Myr, i.e., somewhat younger than the P Pic moving 
group , another nearby young association fe.g.. ITorres et all 
l2008f ). 

Despite being older than most prototypical cTTSs, 
TW Hya is still surrounded by a relatively massive accre- 
tion disc visible at IR and radio wavelengths, and sub- 
ject t o intense scrut iny in the last decade (e.g., iQi et al] 
|2004 iThi et alJfeoiCh . With a typical mass of a few Jupiter 
masses, this disc is found to be inclined at ~10° to the plane 
of the sky; it consists at least of an optically thin inner cav- 
ity extending outwards up to about 4 AUs and of an opti- 
cally thick out er disc extending to distances of up to about 
200 AUs (e.g.. lCalvet et alj|2002t lAkeson et alj|201ll ). sug- 
gesting that planet formation is ongoing. The gas-to-dust 
mass in the disc may be lower than the standard interstellar 
value, suggesting that a significant fraction of the primordial 
disc gas has been evaporat ed due to the s trong X-ray flux 
from the central protostar (|Thi et al.ll2010l ). 

The photospheric temperature of TW Hya estimated 
from th e strength of spectral lines is found to be around 
4150 K i|Torres et al.ll2003l ; lYang et al.ll2005h . only slightly 
hotter than the value of 4000 K t hat the V — R and 
V — I photometric colors indicates (|Rucinski fc Krautterl 
1 19831 ; iRucinski et alj|2008l ); we therefore choose a tempera- 
ture of 4075 ±75 K as a co mpromise. From this temperature, 
evolutionary models (e.g-. ISiess et alJl200Ch indicate that the 
parent protostar (whose contraction is still mostly isother- 
mal at an age of 8-10 Myr) has a typical mass of 0.8 Mq, 
a logarithmic luminosity (with respect to the Sun) of -0.45 
and radius of 1.1 Kq. This agrees well with the brightest 
photometric level re corded for TW Hya (mv = 10.6, e.g., 
IRucinski etail2 008 ) , the distance modulus derived from the 
Hipparcos distance (equal 3.75 ± 0.25) and the bolometric 
correction correspo nding to the phot ospheric temperature 
(-0.95 ± 0.05, e.g.. lBessell et al.lll998r ). yielding a bolomet- 
ric magnitude of 5.90 ± 0.25, i.e., a logarithmic luminosity 
(with respect to the Sun) of -0.46 ± 0.10 (see Fig. [[]). Evo- 
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Figure 1. Observed (open square and error bars) location of 
TW Hya in the HR diagr am. The PMS ev olutionary tracks and 
corresponding isochrones jSiess et al.ll200d) assume solar metal- 
licity and include convective overshooting. 



lutionary models also indicate that, at this age, TW Hya is 
likely no longer fully-convective and hosts a radiativ e core 
extending up to ~0.35 R* in radius (|Siess et al.ll2000l ). 

The spectrum of TW Hya was recently found to un- 
dergo regular radial velocity (RV) fluctuations with a pe- 
riod of 3.56 d, first attributed to a close-in giant planet 
l|Setiawan et al.ll2008l ) then to the presence of dark spots at 
the surface of the star (better reproducing in particular how 
the R V amplitude varies with wavelength, iHuelamo et alj 
2008). It therefore strongly suggests that the period of the 
RV fluctuations is the rotation period of TW Hya, imply- 
ing an equatorial rotation velocity of 15.6 kms -1 (given 
the assumed radius of 1.1 R0, see above). Numerous es- 
timates of the line-of-sight projected equatorial rotation ve- 
locity (noted vsini where i is the inclination of the rota- 
tion axis to the line of sight) are available in the literature, 
ranging from 4 to 15 kms" 1 (|Torres et al.ll2003l ): the most 
reliable of these estimates (those including a basic mod- 
elling of the micro- and/or macro-turbulent broadening of 
the spectral lines in particular) are all on the low side of 

Alencar fc Batalhal 



this range (between 4 and 6 kms , e.g. 



120021 ; iTorres et al.ll2003l ; lYang et al.ll2005l ). consistent with 
the fact that TW Hya is seen mostly pole-on (as the disc 
itself). Our own estimate (4 ± 1 kms -1 , see Sec. O implies 
i ~ 15°, in agreement with previous estimates and w ith the 
inclination of the disc (e.g.. lAlencar fc "B atalha 2 0021 ). 

From both magnetic broadening of spectral lines and 
Zeeman signatures of emission lines used as accretion prox- 
ies, kG magnetic fields were recently re ported at the sur- 
face of TW Hya (|Yang et al.ll2005l . 120071 ) . making it an ob- 
vious target for studies of magnetospheric accretion pro- 
cesses. High-resolution X-ray spectra of TW Hya exhibit 
numerous spectral lines formed at low and high tempera- 
tures and allowing o ne to probe both the accretion shock 
and t he corona (e.g., iKastner et al. 2002; iBrickhouse et al.l 
|2010| ); the low-temperature component of the X-ray spec- 
trum (related to the accretion process) is particularly promi- 
nent on TW Hya (with respect to other cTTSs) making it 
an obvious laboratory for studying the physics of magneto- 



spheric accretion processes. Observations are found to agree 
reaso nably well with the predictions of theoretical models 
(e.g.. lGiinther et al.ll2007l) . 

From the equivalent widths and the corresponding line 
fluxes of the optical emission lines (taken as accretion prox- 
ies) and using the published empirical correla tions between 
lines and accretion fluxes l|Fang et al.l 120091 ). we can de- 
rive an estimate of the logarithmic mass accretion rate (in 
units of Mq yr -1 ) at the surface of TW Hya, that we find 
to be equal to —8.9 ± 0.4 at both epochs (see Sec. [3]), in 
agre ement with indepe ndent estimates from optical prox- 
ies (jCurran et al.l l201ll ). As usual, this is larger than the 
estimate derived f rom X-ray data (equal to —9.74 ± 0.05, 
ICurran et al]|201ll ). 

Optical veiling, i.e., the apparent weakening of the pho- 
tospheric spectrum (presumably caused by accretion) , is sig- 
nificant in TW Hya and typically ranges from 30 to 100% 
during both observi ng runs, in agreem ent with previous pub- 



igr_ 

lished reports (e.g.. lYang et al.ll2005l ) 



3 OBSERVATIONS 

Spectropolarimetric observations of TW Hya were collected 
at two different epochs, first from 2008 March 15 to 28, 
then from 2010 February 23 to March 08, using the high 
resolution spectropolarimeter ESPaDOnS on the Canada- 
France Hawaii Telescope (CFHT). ESPaDOnS collects stel- 
lar spectra spanning the entire optical domain (from 370 to 
1,000 nm) at a resolving power of 65,000 (i.e., 4.6 kms -1 ) 
and with a spectral samp ling of 2.6 km s -1 , in either circu- 
lar or linear polarisation (|Donatill2003l ). In 2008, a total of 
13 circular polarisation spectra were collected over a times- 
pan of 14 nights, but with fairly irregular time sampling; 
in 2010, another set of 13 circular polarisation spectra were 
collected, again over 14 nights, but this time with at a much 
more regular rate (of about 1 spectrum per night). All polar- 
isation spectra (except the first 3 in 2008 March, affected by 
a setup problem) consist of 4 individual subexposures (each 
lasting 815 s and 840 s in 2008 and 2010 respectively) taken 
in different polarimeter configurations to allow the removal 
of all spurious polarisation signatures at first order. 

All raw frames are processed with Libre ESpRIT, 
a fully automatic reduction package/pipeline available at 
CFHT. It automatically performs optimal extraction of ES- 
PaDOnS unpolarized (Stokes /) and circularly polarized 
(S tokes V) spectra gro ssly following the procedure described 
in iDonati et al.l l|l997l ) . The velocity step corresponding to 
CCD pixels is about 2.6 kms -1 ; however, thanks to the fact 
that the spectrograph slit is tilted with respect to the CCD 
lines, spectra corresponding to different CCD columns across 
each order feature a different pixel sampling. Libre ES- 
pRIT uses this opportunity to carry out optimal extraction 
of each spectrum on a sampling grid denser than the origi- 
nal CCD sampling, with a spectral velocity step set to about 
0.7 CCD pixel (i.e. 1.8 kms -1 ). All spectra are automati- 
cally corrected of spectral shifts resulting from instrumen- 
tal effects (e.g., mechanical flexures, temperature or pres- 
sure variations) using telluric lines as a reference. Though 
not perfect, this procedure provides spectra with a relative 
RV pr ecision of better than 0.030 kms -1 (e.g.. IDonati et ail 
20083). 
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Table 1. Journal of observations collected in 2008 March and 
2010 March. Each observation consists of a sequence of 4 subexpo- 
sures (each lasting 815 s and 840 s in 2008 and 2010 respectively), 
except for the first 3 in 2008 March for which only 2 subcxposures 
(of 815 s each) could be used (following an instrument setup prob- 
lem). Columns 1 — 4 respectively list the UT date, the Heliocentric 
Julian Date (HJD) and UT time (both at mid-exposure), and the 
peak signal to noise ratio (per 2.6 kms — 1 velocity bin) of each ob- 
servation. Column 5 lists the rms noise level (relative to the unpo- 
larized continuum level I c and per 1.8 kms -1 velocity bin) in the 
circular polarization profile produced by Least-Squares Decon- 
volution (LSD), while column 6 indicates the orbital/rotational 
cycle associated with each exposure (using the ephemeris given 
by Eq.QJ. 
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The peak signal-to- noise ratios (S/N, per 2.6 kms - 
velocity bin) achieved on the collected spectra range between 
100 and 260 depending on weather/seeing conditions (and 
on the instrumental setup), with a median value of about 
240. The full j ournal of observ a tions is presented in Table [T] 

Following iHuelamo et all (|2008l ). we assume that the 
period of the RV variations is the rotation period. Rota- 
tional cycles E are computed from Heliocentric Julian Dates 
(HJDs) according to th e ephemeris (note the different initial 
HJD than that used by IHuelamo et alJliooi ): 



HJD = 2454540.0 4- 3.5683B. 



(1) 



Coverage of the rotation cycle is only moderate in 2008 
March (with a phase gap of about 1/3 of the cycle) and good 
in 2010 March (see Table [TJ. As a result of the difference 
in the assu med initial HJD, our rotation cycles are ahead 
of those of IHuelamo et all (|2008h by 2.774 cycles (at both 




Velocity (km/s) 

Figure 2. LSD circularly-polarized (Stokes V) and unpolarized 
(Stokes /) profiles of TW Hya (top/red, bottom/blue curves re- 
spectively) collected on 2010 February 26 (cycle 1994-1.101). A 
strong Zeeman signature (with a full amplitude of 1.8%) is de- 
tected in the LSD Stokes V profile, in conjunction with the un- 
polarised line profile. The mean polarization profile is expanded 
by a factor of 10 and shifted upwards by 1.1 for display purposes. 



epochs); maximum RVs should thus occur around phase 0.1 

in our ephemeris. 

Least-Squares Deconvolution (LSD. lDonati et al.lll997l ) 
was applied to all observations. The line list we employed 
for LS D is computed from an Atlas9 LTE model atmo- 
sphere |Kurucj|l993l ) and corresponds to a K7V spectral 
type (Teff = 4,000 K and logg = 4.5) appropriate for 
TW Hya. Only moderate to strong atomic spectral lines 
(with line-to-continuum core depressions larger than 40% 
prior to all non-thermal broadening) are included in this list; 
spectral regions with strong lines mostly formed outside the 
photosphere (e.g., Balmer, He, Call H, K and IRT lines) 
and/or heavily crowded with telluric lines were discarded. 
Altogether, about 6,500 spectral features (with about 40% 
from Fe i) are used in this process. Expressed in units of 
the unpolarized continuum level I c , the average noise lev- 
els of the resulting Stokes V LSD signatures range from 
1.9 to 6.4xl0 -4 per 1.8 kms" 1 velocity bin (median value 
2.4xl0 -4 ). Zeeman signatures are detected at all times in 
LSD profiles (see Fig. [2] for an example) and in most accre- 
tion proxies (see Seed]). 



4 SPECTROSCOPIC VARIABILITY 

In this section, we present a basic description of how the 
spectral lines of TW Hya, and in particular the photospheric 
LSD profiles and the selected accretion proxies (i.e., the 
Ca n infrared triplet or IRT, the He I D3 line and the first 
2 lines of the Balmer series), vary with time and with rota- 
tion phase. The idea is to look at how the equivalent widths, 
the RVs and the average magnetic fluxes of the various pro- 
files are modulated with rotation, to get a rough, intuitive 
understanding of how the large-scale field is structured and 
oriented, and how cool photospheric spots and hot chro- 
mospheric accretion regions are distributed over the stellar 
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Figure 3. Rotational modulation of the RV (top row), veiling (second row) and longitudinal field (bottom row) derived from the LSD 
photospheric profiles of TW Hya in 2008 March (left panels) and in 2010 March (right panels). Data collected during rotational cycles 0, 
1, 2 and 3 are respectively shown with red circles, green squares, dark-blue triangles and light-blue stars. Fits with sine/cosine waves are 
included (and shown as dashed lines) to outline (whenever significant) the amount of variability attributable to rotational modulation. 
±1 a error bars on data points are also shown. 



surface - to be compared with the results of the full tomo- 
graphic imaging analysis carried out in the following section. 



4.1 LSD photospheric profiles 

As shown on Fig. [3] (top and middle panels), the photo- 
spheric Stokes I LSD profiles exhibit significant variabil- 
ity at the two observing epochs, both in position (RV) and 
strength (as a result of veilingjj. Showing no more than a 



moderate level of intrinsic variability, RVs are clearly mod- 
ulated by rotation, varying by ~0.5 kms - peak to peak 
around a mean of 12.55 ± 0.10 kms -1 and reaching a max- 
imum at ph ase ~0.1 (correspond ing to phase 0.88 in the 
ephemeris of lHuelamo et af1l2008l). Thi s is fully compatible 
with the results of Huelamo et al.ll|200Sh . indicating that the 
RV variability we repor t is of the sam e origi n and nature 
than that discussed by iHuelamo et all (|2008l ). i.e., caused 
by a cool photospheric spot on TW Hya located at high lat- 
itudes and at phase ~0.35 (0.25 cycle after RV maximum); 



We assume a conservative error bar of ±0.05 kms - and ±2% 
on the relative RV and the relative strength of photospheric LSD 
and Ca II IRT profiles; in particular, the error bar on the RVs of 
LSD profiles is larger than the formal error bar of 0.03 km s _1 , but 



compa rable to the typical RV dispersion quoted bv lHuelamo et all 
j2008h . 
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2008 March, Ca II IRT 



2010 March, Ca II IRT 
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Figure 4. Same as Fig. [3] for the Ca II IRT LSD profiles. 



moreover, it suggests that this spot is apparently persistent 
on timescales of at least 3 yrs. 

LSD profiles of TW Hya are also significantly shallower 
than those expected from a standard star of similar spectral 
type, indicating that the spectrum is significantly veiled (as 
a likely result of accretion), at levels varying from 30% to 
100% (with an average of about 60% ). This is again ty pical 
to previously published reports fe.g. jYang et al.ll2005t ). We 
find that veiling can be mainly modulated by rotation at 
times, as in 2008 March where maximum veiling is observed 
to occur at phase 0.6 (i.e., 0.25 cycle later than the phase at 
which the cool spot is best visible); at other times, veiling is 
dominated mostly by intrinsic variability, as in 2010 March 
where no clear trend with rotation phase is detected (see 
middle panels of Fig. [3]) . 

As mentioned already, Zeeman signatures are detected 
at all times in Stokes V LSD profiles of TW Hya, show- 
ing a fairly canonical shape (i.e., antisymmetric with re- 



spect to the line centre) and an unusually large amplitude 
of almost 2% peak to peak (see, e.g., Fig. [5]). The line-of- 
sight projected component of the field averaged over the 
visible stellar hemisphere and weighted by brightness inho- 
mogeneities (called the longitudinal field and es timated from 
the fi rst moment of the Stokes V profile, e.g.. iDonati et all 
Il997t) ranges from 380 to 700 G, with typical error bars of 
~15 G. Note that our longitudinal field estimates are sig- 
nificantly larger than those fou nd in previously published 
studies (reporting values <150 G. lYang et al.ll2007h ; this dis- 
crepancy can potentially be attributable to long-term evo- 
lution of the large-scale magnetic topology (as suggested 
for a similar misma tch in the case of the cTTS BP Tau, 
IDonati et alj|2008bl ). as the data of lYang et ah! (120071 ) were 
collected about a decade earlier than ours (in 1999). 

As shown on Fig. [3] (bottom panels), the longitudinal 
field curves at both epochs include a significant amount of 
intrinsic variability (especially in 2010 March), but never- 
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theless suggest a ±10-20% level of rotational modulation 
about the mean value (of ~500 G) with maximum longi- 
tudinal field occurring at phase ~0.35, i.e., when the cool 
spot is best visible. Most likely, this intrinsic variability of 
the longitudinal field does not reflect a true change in the 
large-scale magnetic topology, but rather relates to modifi- 
cations in the line formation region and related properties 
(e.g., veiling, cool spot location and size). 

4.2 Ca n IRT emission 

As in previous papers, we use core emission in Ca II IRT 
lines as our main proxy of surface accretion even though 
it usually features a significant contribution f rom the non- 
accre t ing ch romosphere; previous work (e.g., iDonati et al.l 
l2010l . 1201 lab demonstrated in particular that rotational 
modulation in Stokes / and V Ca 11 IRT profiles of cTTSs 
correlates well with that derived from the more conventional 
He 1 D3 accretion proxy (see below), even at fairly low mass- 
accretion rates. To recover the emission component of the 
Ca 11 IRT lines, we start by constructing a LSD-like weighted 
average of the 3 IRT lines; we then subtract the underlying 
(much wider) photospheric absorption profiles, with a single 
Lorentzian fit to the far line wings f see IDonati et al1l2011bl . 
for an illustration). 

The average Ca 11 emission profiles are centred at a RV 
of ~13.40 kms" 1 , i.e., slightly red-shifted (by 0.85 kms" 1 ) 
with respect to the photospheric LSD profiles (see Fig. U 
top panels); this is usual for cTTSs and further confirms 
that Ca 11 emission probes slowly-moving plasma in the post- 
shock regions of accretion funnels, and located close to the 
surface of the star. We also observe that the RVs of the Ca 11 
emission profiles show, in addition to a weak level of intrinsic 
variability, a small amount of rotational modulation, with an 
amplitude comparable to (in 2010 March) or smaller than 
(in 2008 March) that seen in the LSD photospheric profiles 
and much smaller (by at least a factor of 15) than the ro- 
tational broadening of spectral lines (vsini ~ 4 kms -1 , see 
Sec. [5J. It suggests in particular the presence of a localised 
region of accretion-powered emission near the pole, and best 
viewed (at both epochs) at a rotational phase of ~0.5 (i.e., 
at mid transit between minimum and maximum RVs). This 
approximately coincides with the phase of maximum veiling 
in 2008 March (~0.6, see above) and with that at which 
the cool photospheric spot (probed through LSD profiles) is 
best visible (~0.35), the associated phase shifts (of +0.10 
and -0.15 rotation cycle respectively) being not necessarily 
very significant at high latitudes and with such a high level 
of intrinsic variability. 

Emission strength is also varying with time (see middle 
panels of Fig. 2}, with equivalent widths typically within the 
range 20-30 kms" 1 (0.055-0.085 nm) and featuring at times 
short-lived sporadic bursts during which emission is tem- 
porarily enhanced by as much as 50% (e.g., at cycle 1.654 
and 1.668 in 2008 March). Rotational modulation is mod- 
erate, potentially up to 30% peak-to-peak in 2008 March 
with maximum emission occurring near phase 0.6 (in con- 
junction with the phase of maximum veiling, see above); in 
2010 March however, rotational modulation is much lower 
(lower than 10% peak-to-peak) and mostly hidden by intrin- 
sic variability. 

As for LSD profiles, clear Zeeman signatures are de- 



tected at all times in conjunction with the emission compo- 
nent of Ca 11 IRT lines, with peak-to-peak amplitudes reach- 
ing up to 10% (in 2010 March). The corresponding longi- 
tudinal fields range between -380 and -730 G (with typical 
error bars of about 30 G), except during the previously men- 
tioned emission burst (when the longitudinal field is reduced 
to —150 G). Note in particular the striking polarity con- 
trast between the (constantly negative) longitudinal fields 
measured from the Ca 11 emission lines and those (always 
positive) derived from the LSD profiles; this clearly demon- 
strates that LSD profiles and IRT lines probe different ar- 
eas (of different magnetic polarities) over the stellar surface 
of TW Hya (as on V2129 Oph and BP Tau, IDonati et ail 
l2011al . l2008bl V The field values we derive from the Ca 11 IRT 
profiles are comparable (both in strength and s ign) to those 
reported in a previous paper (|Yang et al.ll2007l ). in contrast 
to longitudinal fields from LSD photospheric profiles (found 
to be significantly stronger in our study, see above); it may 
suggest that the potential field evolution between the two 
studies concerns the magnetic topology as a whole and not 
just its overall strength. 

Again, intrinsic variability is dominant at both epochs 
and no more than a moderate level of rotational modulation 
is observed, with strongest longitudinal fields occurring ap- 
proximately when the cool photospheric spot is best visible 
(i.e., around phase 0.35); as for longitudinal fields estimated 
from LSD profiles, we suggest that this intrinsic variability 
reflects changes in the line formation region (and in partic- 
ular in the accretion region) rather than in the large-scale 
magnetic topology itself. 

4.3 He 1 D3 emission 

The He I D3 line, used in most cases as the reference accre- 
tion proxy, exhibits, in the spectrum of TW Hya, an emission 
profile that is strongly asymmetric with respect to the line 
centre (see Fig. [SJ. This profile can be decomposed into a 
narrow central emission centred at ~18.2 kms -1 (i.e., red- 
shifted by about 5.6 kms -1 with respect to the LSD profiles) 
and whose full-width at half maximum is ~36 kms" 1 , and 
a broad emission component (full width at half maximum 
~110 kms -1 ) red-shifted by typically 20-30 kms -1 with re- 
spect to the narrow component; the equivalent width of the 
narrow component (35 and 45 kms -1 on average in 2008 
March and 2010 March respectively, or equivalently 0.07 
and 0.09 nm) is about half that of the broad component 
(~75 kms -1 or 0.15 nm). Zeeman signatures with typical 
amplitudes of ~20% are observed in conjunction with the 
narrow emission component, but not with the broad com- 
ponent (see Fig. [5]). Their shape strongly departs from the 
usual antisymmetric pattern (with respect to the line cen- 
tre) likely indicating that the li ne forms in a region with 
significant velocity gradients fe.g. JSanchez Almeida fc Litesl 
1 1992h . We suggest that this narrow emission component 
comes from the postshock region in which the postshock 
plasma is strongly decelerating though still falling towards 
the stellar surface (hence the red-shifted emission and the 
non-antisymmetric Zeeman signature). Similarly, we pro- 
pose that the broad emission component probes the whole 
accretion funnels rather than just the postshock region 
(hence the larger line width and red-shift, reflecting the 
higher plasma velocities involved). In the following, we focus 
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Figure 5. Variations of the unpolarized (Stokes /, left panel) and circularly-polarized (Stokes V, right panel) profiles of the He I Dz 
emission of TW Hya in 2008 March (left columns of both panels) 2010 March (right columns). The unpolarized emission profile is 
asymmetric with respect to the line centre and can be described as the sum of a central narrow component and a broad (often red- 
shifted) component (see text). A strong Zeeman signature (with a full amplitude of ~20%) is detected in conjunction with the central 
narrow He I emission, but not with the broad red-shifted component. To emphasize variability, the average profile over each run is shown 
in red. Rotation cycles (as listed in Table 1) and 3a error bars (for Stokes V profiles only) are shown next to each profile. 



on the narrow emission component mostly, being the most 
relevant one for our study on the large-scale field topology, 
photospheric spots and accretion regions at the surface of 
TW Hya. 

The temporal variability of the narrow He I emission is 
shown in Fig. |S]for both observing epochfl We find that the 
RV of the narrow He I emission is modulated by rotation in 
a similar way than the Ca II IRT emission, suggesting again 
the presence of an accretion region located around phase 
0.3-0.4 (at mid distance between minimum and maximum 
RVs), i.e., more or less cospatial with the cool photospheric 
spot probed through LSD profiles (see above). The semi- 
amplitude of the RV variations reaches 0.5-0.8 kms -1 for 
the narrow He I emission (see Fig. [6] top panels) , signifi- 
cantly larger than those seen in the Ca II IRT lines; this is 
in agreement with the much smaller dilution that the He I 
line suffers from the non-accreting chromosphere. Assuming 
that this dilution is negligible for the He I line and thus that 



2 We assume an error bar of ±0.25 kms" 1 and ±2.5 kms" 1 
on the relative RV and the relative strength of the narrow He I 
emission component, i.e., 5x larger than those used for the Call 
IRT profiles to reflect the reduced accuracy resulting from the 
subtraction of the broad He I emission component. 



the associated RV fluctuations directly track the Doppler 
motion of the hot spot, we infer that this accretion region 
is off-centred with respect to the pole by typically 10° (for 
wsini=4 kms -1 ). 

We also find that the equivalent width of the narrow 
He I emission is varying with time. Rotational modulation is 
again moderate, up to about 30% and peaking at phase ~0.4 
in 2008 March and no larger than 10% in 2010 March (see 
Fig. [6] middle panels); intrinsic variability is significant and 
dominates in 2010 March. These conclusions mostly confirm 
what the Ca II IRT lines already suggested (see above) . Sur- 
prisingly, the emission burst detected in Ca II does not seem 
to show up in the narrow He I emission component, although 
clearly detected in the He I line (broad component); this may 
however reflect the limits of our two-component decompo- 
sition (used to describe the He I line) in the particular case 
of emission bursts (where both components happen to have 
roughly equal RVs, see Fig. [5j . 

The longitudinal fields we derive from the observed Zee- 
man signatures range from -2 to -3.5 kG (see Fig. HJ] bot- 
tom panels) when we take into account that the narrow He I 
emission is the only contributor. Unsurprisingly, the field po- 
larity is the same as that derived from Ca II IRT emission, 
demonstrating that both lines are reliable probes of accre- 
tion regions (even though Ca II emission suffers a strong di- 
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Figure 6. Same as Fig. [3] for the narrow component of the He I D3 profiles. 



lution from the non-accreting chromosphere) . Our He I lon- 
gi tudinal field estim ates are significantly larger than those 
of lYang et alj (|2007) by typically a factor of 2; this discrep- 
ancy is however only apparent and reflects their assumption 
that the full He I emission (including the broad emission, 
rather than the narrow component only) contributes to the 
Zeeman signatures. 



As for the Ca n IRT emission, intrinsic variability domi- 
nates and only moderate rotational modulation is observed, 
with maximum longitudinal field (of about -3 kG) occur- 
ring approximately when the cool photospheric spot is best 
visible; once again, we suspect that frequent changes in the 
properties of the accretion region are causing the enhanced 
dispersion of the observed longitudinal fields. 



4.4 Balmer emission 

Balmer lines are also strongly in emission in the spectrum 
of TW Hya, with Ha and H/3 reaching typical equiva- 
lent widths of 8,000-10,000 kms" 1 (18-22 nm) and 1,700- 
2,100 kms -1 (2.8-3.4 nm) respectively. At times, they show 
a conspicuous blue-shifted absorption component (e.g., at 
cycles 0.274-0.546 in 2008 March) indicating the presence 
of a strong wind escaping the star at an average projected 
velocity of ~50 kms - . 

Both lines exhibit significant fluctuations with time (see 
Fig. [7] however, since the profiles and emission strengths of 
Ha and H/3 do not convincingly repeat between consecutive 
rotation cycles (e.g., cycles 0.517 and 2.515 in 2008 March, 
cycles [199+] 0.255 and 2.221 in 2010 March), these vari- 
ations cannot be attributed to rotational modulation and 
are more likely caused by intrinsic temporal evolution of the 
accretion and wind flows. In 2008 March, the blue-shifted 
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wind absorption is causing most of the observed variability 
in Ha and H/3; in 2010 March, the emission fluxes of both 
Balmer lines correlate well with that of the broad (and not 
of the narrow) He I emission component. 

We also note that the red wing of H/3 is depressed 
at times with respect to the average profile (e.g., at cycle 
[199+] 0.805 in 2010 March). However, these depressions 
never show up as true absorption as they r egularly do on 
some cTTSs (e.g., AA Tau and V2129 Oph, iBouvier et all 
l2007bl ; iDonati et all 1201 lal l: we thus think it is ambiguous 
and potentially misleading to relate these events to the cross- 
ing of accretion funnels in front of the visible hemisphere of 
TW Hya (as done for AA Tau and V2129 Oph), as they 
could just as well result from varying projection effects of 
the magnetospheric funnel structure. 



— 2.1, —1. 6 and —1.2 respe ctively (using empirical correla- 
tions from lFang et al . 2009). We thus conclude that the aver- 
age logarithmic mass accretion rate of TW Hya (in Mq yr _1 ) 
is equal to —8.9 ± 0.4 at both epochs, in rough agreement 
with the estimate of lCurran et all (|201ll ). equal to -9.2T0.2 
and also derived from optical lines (though from fully inde- 
pendent data and a different analysis). 

Mass accretion rates can in principle also be estimated 
(though less acc urately) through the full width of Ha at 
10% height (e.g- lNatta et alfeOO^ICieza et al.ll2(5lol ) . In the 
case of TW Hya, Ha shows a full width of ~400 kms -1 on 
average, implying a logarithmic mass accretion rate estimate 
of —9.0 ± 0.6 (in Mq yr _1 ) in agreement with our previous 
estimate. 



4.5 Mass-accretion rate 

From the average equivalents widths of the Ca II IRT, He I 
(narrow component) and H/3 lines, we can derive logarith- 
mic line fluxes (with respect to the luminosity of the Sun 
Lq), equal to —4.8, —4.8 and —3.3 respectivehfl implying 
logarithmic accretion luminosities (with respect to Lq) of 

3 To derive line fluxes from normalized equivalent widths, we 
approximate the continuum level by a Planck function at the 
temperature of the stellar photosphere. 



5 MAGNETIC MODELLING 
5.1 Overview of the method 

Now that the temporal variability of photospheric LSD and 
Ca II IRT emission profiles has been examined in detail, we 
aim at converting the corresponding data sets into maps of 
the large-scale magnetic topology, as well as distributions of 
surface cool spots and of chromospheric accretion regions, at 
the surface of TW Hya. To achieve this, we apply our tomo- 
graphic imaging tech nique, described extensive ly in previous 
similar studies (e.g., IDonati et al]|2010l . l2011al ). 
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Following the principles of maximum entropy, our code 
automatically and simultaneously derives the simplest mag- 
netic topology, photospheric brightness image and accretion- 
powered Ca II emission map compatible with the series of 
observed Stokes / and V LSD and Ca II IRT profiles. The 
reconstruction process is iterative and proceeds by compar- 
ing at each step the synthetic Stokes I and V profiles corre- 
sponding to the current images with those of the observed 
data set. The magnetic topology is described through its 
poloidal and toroidal component s expressed as sph erical- 
harmonics (SH) expansions (e.g.. iDonati et ail 120061. The 
spatial distributions of photospheric brightness (with re- 
spect to the quiet photosphere) and those of accretion- 
powered Ca II emission (in excess of and with respect to 
that produced by the quiet chromosphere) are modelled as 
series of independent pixels (typically a few thousand) on a 
grid covering the visible surface of the star (with spots in 
the brightness image assumed to be darker/cooler than the 
quiet photosphere and accreting regions in the Ca II emission 
maps supposed to be brighter than the quiet chromosphere). 

Synthetic profiles are computed by summing up the 
elementary spectral contributions from all image pixels 
over the visible stellar hemisphere, taking into account the 
relevant local parameters of the corresponding grid cells 
(e.g., brightness, accretion-powered excess emission, mag- 
netic field strength and orientation, radial velocity, limb an- 
gle, projected area). Since the problem is partly ill-posed, 
we stabilise the inversion process by using an entropy cri- 
terion (applied to the SH coefficients and to the bright- 
ness/excess emission image pixels) aimed at selecting the 
magnetic topology and images with minimum information 
among all those compatible with the data. The relative 
weights attributed to the various SH modes can be imposed, 
e.g., for purposely producing antisymmetric or symmetric 
field topologies with respect to the centre of the star (b y 
favouring odd or even SH modes. IDonati et al 1 l2007l . l2008bll . 
More details concerning the specific de scription of lo c al pro - 
files used in the model can be found in IDonati etafl (|20ld l. 

5.2 Application to TW Hya 

Our imaging model assumes that the observed profile vari- 
ations are mainly due to rotational modulation; all other 
sources of profile variability (and in particular intrinsic vari- 
ability) cannot be properly reproduced and thus contribute 
as noise into the modelling process, degrading the imaging 
performance and potentially even drowning all relevant in- 
formation. Filtering out significant intrinsic variability from 
the observed profiles is thus worthwhile to optimise the be- 
haviour and convergence of the imaging code. 

We implement this by applying specific corrections on 
our data set. We first suppress veiling by scaling all LSD 
Stokes / and V photospheric profiles, to ensure that un- 
polarized lines have the same equivalent widths. We also 
remove the non rotationally-modulated part in the observed 
fluctuations of Ca II IRT emission, by fitting them with 
a sine+cosine wave (see Fig. [5| and by scaling the corre- 
sponding Stokes / and V profiles accordingly, thus ensur- 
ing that equivalent widths of unpolarized lines match the 
optimal fit. Although only approximate (especially the re- 
moval of the intrinsic variability), this procedure is at least 
straightforward and has proved successful when applied to 



real data (e.g.. IDonati et aill2011al l and efficient at retain- 
ing rotational modulation mostly. In the particular case of 
the present data set, the most significant effect of this filter- 
ing process is to remove the emission burst detected in 2008 
March (at cycle 1.654 and 1.668) by scaling down the 2 cor- 
responding Ca II IRT profiles (shown as green open squares 
on Fig. [8] left panel); the remaining modulation left in the 
Ca II profiles are the ~30% and < 10% fluctuations reported 
in Sec. H 

We use Unno-Rachkovsky's equations known to pro- 
vide a good description of the local Stokes I and V 
profiles (including magneto-optical effects) in the pres- 
ence of both weak and strong m agnetic fields (e.g., 
lLandi degl'Innocenti fc Landomll2004l , Sec. 9.8) despite their 
being based on the assumption of a simple Milne-Eddington 
atmospheric model. The model parameters used in the spe- 
cific case of TW Hya are mostly identical to th ose used in 
our previous studies (IDonati et ail |2010| . 1201 fal l. The emis- 
sion profile scaling factor e, describing the emission enhance- 
ment of accretion regions over the quiet chromosphere, is 
once again set to e = 10. We also assume (as in previous 
studies) that the magnetic topology of TW Hya is antisym- 
metric with respect to the centre of the star; such field con- 
figurations are indeed best suited to explain ho w accretion 
regio ns can form at high latitudes mostly (e.g., lLong et ail 
2008), as often the case for cTTSs. This assumption has 
little impact on the reconstructed magnetic maps over the 
visible regions of the stellar surface and allows a consistent 
and realistic (though not necessarily exact) field distribution 
to be recovered on the unseen hemisphere as well. 

The magnetic, brightness and accretion maps we recon- 
struct for TW Hya at both epochs are shown in Fig. with 
the corresponding fits to the data shown in Fig. 1101 The 
SH expansion describing the field was limited to terms with 
£ < 5; attempts with I = 10 indicate that little power is re- 
constructed in higher order modes (£ > 6), reflecting essen- 
tially the limited spatial information accessible to Doppler 
tomography at low v sin i values. Error bars on Zeeman 
signatures were artificially expanded by a factor of 3 (for 
LSD profiles) and 2 (for Ca II emission) at both epochs to 
take into account the significant level of intrinsic variability; 
moreover, Zeeman signatures corresponding to the 2 Ca II 
profiles recorded during the emission burst of 2008 March 
(at cycles 1.654 and 1.668) were excluded from the fit, being 
strongly discrepant with others taken at very similar phases 
(but different cycles, e.g., 3.633). The fits we finally obtain 
correspond to a reduced chi-square Xi of 1.2, starting from 
initial values of 59 and 47 for the 2008 March and 2010 
March data sets respectively (with scaled-up error bars on 
Zeeman signatures). 

As a by-product, we obtain new estimates for various 
spectral characteristics of TW Hya; in particular, we find 
that the RV is equal to 12.55 ± 0.10 kms -1 , that the mean 
red-shift of Ca II emission profiles with respect to LSD pro- 
files is 0.85 ± 0.1 kms -1 , and that the v sin i is 4± 1 kms" 1 . 
Finally, we find that the profiles are best fitted for a value of 
the local filling factor ip, describing the relative proportion 
of magnetic areas at any given point of the stellar surface, 
equal to ~0.4. 
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Figure 8. Measured (open symbols) and fitted (pluses) equivalent widths of the Ca II IRT LSD profiles of TW Hya in 2008 March (left 
panel) and 2010 March (right panel). The model wave (dashed line) providing the best (sine+cosine) fit to the data presumably traces 
rotational modulation, while the deviation from the fit illustrates the strength of intrinsic variability. The open symbols are defined as 
described in Fig. [3] 



5.3 Modelling results 

The reconstructed large-scale magnetic topology of TW Hya 
is very similar at both epochs, featuring a strong and mostly 
axisymmetric field as expected from the small amount of ro- 
tational modulation that Zeeman signatures exhibit. The 
magnetic field at the surface of the star has an average 
intensity of ~1.5 kG, concentrating mostly in the radial 
and meridional components; in particular, the radial field 
map consists of a spot of negative radial field near the pole 
(where the field reaches up to 3 kG, in agreement with 
the field strengths estimated from the narrow He I emis- 
sion component) surrounded with a ring of positive radial 
field at low latitudes (see Fig. |9}. We find that the field 
is almost purely poloidal and mostly axisymmetric at both 
epochs, with < 5% of the reconstructed magnetic energy 
in the toroidal and non-axisymmetric (i.e., SH modes with 
m > 1/2) poloidal fields. 

The octupolar component of the poloidal field largely 
dominates at both epochs (reaching polar strengths of 2.5 
and 2.8 kG in 2008 March and 2010 March respectively) as 
expected from the topologies of the reconstructed radial and 
meridional field maps; the dipole component is typically 4-6 
times weaker, with a polar strength of 0.4 and 0.7 kG in 2008 
March and 2010 March respectively. The octupole is aligned 
with respect to the rotation axis to better than 10°, with 
the negative pole being visible from the Earth. The recon- 
structed dipole component is also aligned with the rotation 
axis in 2010 March but tilted by about 45° (towards phase 
0.5) in 2008 March. However, given the relative weakness 
of the dipole field and the incomplete phase coverage of the 
dataset at this epoch, it is unclear whether this tilt is signif- 
icant; extensive simulations indeed suggest that retrieving a 
weak dipole component in a predominantly octupolar large- 
scale field is uncertain, especially at low usini's. At both 
epochs, the visible pole of the dipolar component is the pos- 
itive pole, suggesting that the dipole and octupole compo- 
nents are essentially anti-parallel in 2010 March. Given the 
high-level of intrinsic variability in the data, we caution that 
the (limited) change in the large-scale magnetic topology of 



TW Hya (especially regarding the dipole component) be- 
tween the two epochs may potentially be spurious. 

The brightness distributions at the surface of TW Hya 
that we recover show a cool spot near the pole at both 
epochs (see Fig. [9]). This spot, covering about 2% of the 
stellar surface, is shifted from the pole by about 10° to- 
wards phase 0.35-0.40, as expected from the RV variations 
of photospheric LSD profiles (see Sec. [4]). This spot is lo- 
cated right on the region of strong negative radial field near 
the pole (the magnetic pole of the octupole component), ex- 
plaining a posteriori why photospheric Zeeman signatures 
are essentially insensitive to this strong negative radial field 
region and probe instead fields of opposite polarities (pos- 
itive radial and meridional fields at intermediate latitudes, 
see Fig. [§]); being very cool, this spot emits much less pho- 
tons (per surface area) than the surrounding photosphere 
and therefore much less polarised photons as well, making 
it essentially invisible to Zeeman signatures of photospheric 
LSD profiles (at least at visible wavelengths). 

The distributions of accretion-powered excess emission 
also show at both epochs a clear spot near the pole, and cov- 
ering 2-3% of the stellar surface (see Fig.[5]). This accretion 
region approximately overlaps with the cool photospheric 
spot, with possibly a small phase shift (towards phase 0.50 
in 2010 March in particular) suggesting that the accretion 
region may slightly trail the photospheric spot. These dis- 
tributions explain a posteriori why Zeeman signatures of 
Ca 11 (and He 1) emission profiles mostly probe fields with 
opposite polarities than those traced by LSD photospheric 
profiles; while Ca 11 (and to a much larger extent He 1) emis- 
sion profiles are mostly sensitive to the polar regions (host- 
ing strong negative radial fields), LSD photospheric profiles 
provide a complementary view of the visible stellar surface 
(and in particular of intermediate latitude regions, hosting 
positive radial and meridional fields). 

The near-polar accretion region apparently also extends 
to low latitudes at phase ~0.6 in 2008 March, potentially in- 
dicating that sporadic non-polar accretion is also occurring 
at times on TW Hya. This low-latitude appendage directly 
reflects the ~30% full-amplitude modulation observed in the 
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Figure 9. Maps of the radial, azimuthal and meridional components of the magnetic field B (first and third rows, left to right panels 
respectively), photospheric brightness and excess Ca II IRT emission (second and fourth rows, first and second panels respectively) at 
the surface of TW Hya, in 2008 March (top two rows) and 2010 March (bottom two rows). Magnetic fluxes are labelled in G; local 
photospheric brightness (normalized to that of the quiet photosphere) varies from 1 (no spot) to (no light); local excess Call emission 
varies from (no excess emission) to 1 (excess emission covering 100% of the local grid cell, assuming an intrinsic excess emission of 
lOx the quiet chromospheric emission). In all panels, the star is shown in flattened polar projection down to latitudes of —30°, with the 
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Figure 10. Maximum-entropy fit (thin red line) to the observed (thick black line) Stokes / and Stokes V LSD photosphcric profiles (top 
panels) and Ca II IRT profiles (bottom panels) of TW Hya. In each panel, the left and right columns correspond to the 2008 March and 
2010 March data respectively. The light-blue curve in the bottom left panel shows the (constant) contribution of the quiet chromosphere 
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emission flux of Ca II (and of He i) profiles peaking around 
phase 0.6 (see Sec.[4| and is consistent with veiling estimates 
(also peaking around phase 0.6 at this epoch). Independent 
confirmation that this non-polar accretion event is real - 
rather than caused by over-interpreting intrinsic variability 
into rotational modulation - comes from the weird Zeeman 
signatures of Ca n (and He i) lines collected during the emis- 
sion burst of 2008 March (at cycles 1.654 and 1.668) and 
excluded from the fitted data set (given their very specific 
shapes, see Fig. I10[) . By subtracting the expected Zeeman 
signatures (the red lines on Fig. I10|) from the observed ones 
(the black lines) for these two profiles, we end up with the 
Zeeman signatures associated with the excess emission spe- 
cific to the emission burst and find that they trace a radial 
field of opposite (i.e., positive) polarity than the main Zee- 
man signatures collected at similar phases but outside of 
the accretion burst (e.g., at cycle 3.633). This unambigu- 
ously demonstrates that sporadic accretion is truly occur- 
ring on TW Hya in regions of positive radial field, i.e., at 
much lower latitudes than where most of the disc material is 
accreted (on the pole), therefore validating the reality of the 
low-latitude appendage to the near-polar accretion region in 
2008 March. 



6 SUMMARY & DISCUSSION 

We present in this paper new results on magnetospheric ac- 
cretion processes taking place at the surfaces of forming Sun- 
like stars; towards this aim, we use spectropolarimetric data 
collected on the evolved cTTS TW Hya at 2 different epochs 
(2008 March and 2010 March) and with ESPaDOnSOCFHT 
in the framework of the MaPP program. Strong Zeeman sig- 
natures are detected in the LSD profiles of photospheric lines 
of TW Hya, as well as from the Ca II and He I emission lines 
commonly used as accretion proxies, revealing the presence 
of longitudinal fields ranging from -3.5 to 0.7 kG depending 
on the set of lines and on the observing epoch; in particu- 
lar, longitudinal fields of opposite polarities are detected (at 
both epochs) from the two sets of lines (i.e., LSD profiles 
and accretion proxies), clearly indicating that they probe 
different regions of the stellar surface. Using our dedicated 
tomographic imaging to ol (tested and validated w ith sev- 
eral similar studies, e.g.. iDonati et all |2010| . l2011al . among 
the most recent ones), we convert these data sets into images 
of the large-scale magnetic topology, as well as distributions 
of the photospheric brightness and of the accretion-powered 
excess emission, at the surface of TW Hya. 

We find that the large-scale magnetic field of TW Hya 
is strong (with a typical average surface intensity of 1.5 kG) 
and mostly poloidal and axisymmetric (with respect to the 
rotation axis). More specifically, the poloidal field is domi- 
nated by a strong octupole component (with a polar strength 
of 2.5-2.8 kG) mostly aligned with the rotation axis (to bet- 
ter than 10°) and exhibiting its negative pole to an Earth- 
based observer. The dipole component is much weaker, by 
typically a factor of 4-6 (with a polar strength ranging from 
~400 G in 2008 March to ~700 G in 2010 March), and is 
found to be mostly anti-parallel with respect to the octupole 
(though apparently tilted at about 45° with respect to the 
rotation axis in 2008 March). We however caution that the 
latter results need further confirmation from additional data 



given the fairly strong octupole to dipole intensity contrast, 
and the significant level of intrinsic variability in the data. 

Regarding its large-scale magnetic topology, TW Hya 
is fairly similar to V2129 Oph, also hosting a strong and 
mostly poloidal, axisymmetric large-scale fi eld in which 
the oc tupole component largely dominates (|Donati et al.l 
l2011al . though with parallel, rather than antiparallel, dipole 
and octupole moments); it is however significantly different 
from AA Tau and BP Tau on the one hand (also hosting 
a strong, mostly axisymmetric, large-scale magnetic field, 
but wit h a larger and som etimes dominant dipole com- 
ponent, IDonati et all |2010T ) and V4046 Sgr, CV Cha or 
CR Cha on the other hand (exhibiti ng a much weaker , 
mostly non-axisymm etric poloidal field. iHussain et aL 2009; 
IDonati et all 1201 lbT ). This is in good agreement with the 
overall trends on the properties of dynamo fields as derived 
from previo us studies, both on cTTSs and main sequence 
dwarfs (e.g.. | Morin et aT. | l2008l ; IDonati fc Landstreetfl2009l : 
M orm et alJbOlOh ;' more specifically, these trends show that 
cool stars (i) are capable of producing intense axisymmetric 
large-scale fields with strong dipole components when fully 
convective (as for A A Tau and BP Tau), (ii) start losing 
their large-scale dipole while keeping a strong axisymmet- 
ric octupole when they start developing a radiative core (as 
for V2129 Oph and TW Hya) and (iii) finally end up with 
weak, mostly non-axisymmetric poloidal fields when their 
radiative cores occupy a large-enough proportion of the star 
(with an outer radius of at least 0.4-0.5 R+, as for V4046 Sgr, 
CV Cha and CR Cha). 

We also confirm that TW Hya hosts a co ol photospheric 
spot n ear the pole, as initially proposed bv lHuelamo et all 
(2008) to account for the observed regular RV modulation 
of its spectrum; we find that this spot is located near the 
pole, shifted (by about 10°) towards phase 0.35-0.40 (in the 
ephemeris of Eq.[T]) and apparently very stable on timescales 
of several years. This spot mostly overlaps with the main 
magnetic pole, which thus becomes essentially invisible to 
LSD photospheric profiles. Most likely, this near-polar, long- 
lived cool spot is a direct consequence of the strong large- 
scale magnetic topology, with the intense fields (of ~3 kG) 
present at the magnetic pole presumably inhibiting outward 
energy transport by convection more than everywhere else 
at the surface of the star. In this respect as well, TW Hya 
resembles V2129 Oph whic h also hosts a similar cool photo- 
spheric spot near the pole (|Donati et al.ll2011al ). 

We finally obtain that TW Hya hosts a near-polar re- 
gion of accretion-powered excess Ca n (and He i) emission, 
mostly overlapping with the cool photospheric spot and co- 
inciding with the main magnetic pole. This confirms in par- 
ticular that accretion occurs mostly poleward at the surface 
of TW Hya, again sim ilar to what is reported for V2129 Oph 
(|Donati et al.j|2011al ). In addition to this regular poleward 
accretion, we also report the detection of sporadic accre- 
tion events towards lower latitudes. These episodic events 
can show up as enhanced rotational modulation in Ca n 
and He I emission (as in our 2008 March images); they can 
even manifest themselves at times as short-lived flaring-like 
emission bursts (as in cycles 1.654 and 1.668) producing 
specific Zeeman signatures of opposite polarities than (and 
thus partially cancelling out with) those associated with the 
main polar accretion region (octupoles exhibiting positive 
and negative radial fields at low and high latitudes respec- 
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tively, or vice versa, within a given hemisphere). We note 
that the epoch at which equator- ward accretion occurs (2008 
March) is also that at which the large-scale dipole field is 
weakest and potentially inclined with respect to the rotation 
axis (towards the phase of maximum Ca II emission approx- 
imately); although this needs confirmation from new data 
with better phase sampling than our 2008 March data set, 
it is at least consisten t with theoretical expectations (e.g., 
iRomanova et al1l201ll ). 

From the emission fluxes observed in spectral lines com- 
monly used as accretion proxies, we conclude that the log- 
arithmic mass accretion rate at the surface of TW Hya is 
equal to —8.9 ± 0.4 (in M©yr _1 ) at both epochs, in rea- 
sonable agreement with r ecent independent determinations 
from optical spectr a fe.g.JCurran et aDl201lh . Following the 
theoretical work of lBessolaz et al.l (|2008T ). we can infer that 
TW Hya is capable of magnetically disrupting its disc up to 
a distance r mag of about 4 ± 1 i?* if the large-scale dipole is 
as strong as 700 G, but only 3 ± 1 JZ* if the dipole strength 
is only 400 G. In the latter case (corresponding to our 2008 
March observations) , this is just barely beyond the distance 
at which the dipole field starts to dominate over the octupole 
field (ranging from 1.7 to 2.2 R+ for octupole to dipo le po- 
lar strength ratios of 4-6, e.g., iGregorv et al-l feoiO): this 
potentially explains why equator-ward accretion events are 
observed preferentially at this epoch. When compared to 
the corotation radius r COT ~ 8.3 ii* (at which the Keple- 
rian period equals the stellar rotation period), we find that 
fmag/fcor potentially ranges from 0.35 ± 0.10 to 0.50 ± 0.10 
depending on the dipole field strength, far below the value 
(of ~ 1) at which star/disc magnetic coupling can start in- 
ducing a si gnificant outflow throu gh a propeller-like mecha- 
nism (e.g.. lRomanova et al.ll2004h and thus force the star to 
spin down and rotate as slowly as one cycle per week. 

We thus speculate that TW Hya, with its relatively 
short rotation period (of ~3.56 d), is in a phase of rapid 
spin-up, mostly reflecting the fact that its large-scale dipole 
field is no longer strong enough to disrupt the disc beyond 
0.5 r cor and fuel the associated leak of angular momentum 
(e.g., through a propeller-like mechanism); this situation is 
similar, though admittedly not exactly ide ntical, to that de- 
scribe d in the recent MHD simulations of lZanni fe Ferreiral 
(2009). In this respect, TW Hya is different from V2129 Oph 
whose rotation period (of ~6.53 d) is still rather long. We 
propose that TW Hya and V2129 Oph are both undergo- 
ing a process of rapid spin-up; however, while V2129 Oph 
is only at the beginning of this phase, TW Hya is in a 
much more advanced state (suggesting at the same time that 
TW Hya has been partly convective for a longer timespan 
than V2129 Oph). 

As MaPP observations accumulate and are being anal- 
ysed, we are progressively building a consistent description 
of what dynamo magnetic fields of cTTSs look like and how 
they impact the formation process, in particular regarding 
the way magnetospheric accretion processes can succeed in 
slowing down the rotation of cTTSs. More observations are 
obviously needed to validate more firmly and extend fur- 
ther the initial conclusions drawn from the existing MaPP 
studies including this one. MaPP is still ongoing, with time 
granted on ESPaDOnSQCFHT (and NARVALOTBL) up to 
semester 2012b; potential extensions, aimed at carrying out 



a longer survey for a few specific cTTSs and/or at expanding 
the sample to wTTSs, are under study at the moment. 
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